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Abstract 

We explore aspects of the phase structure of SU (2) and SU (3) lattice gauge 
theories at strong coupling with many flavours Nf of Wilson fermions in the fun- 
damental representation. The pseudoscalar meson mass as a function of hopping 
parameter is observed to deviate from the expected analytic dependence, at least 
for sufficiently large Nf. Implications of this effect are discussed, including the 
relevance to recent searches for an infrared fixed point. 



PACS numbers: 11.5.15.Ha 



1. Introduction 



There has been a lot of interest recently in identifying an infrared fixed point, and the 
possibility of an associated conformal window, for SU(N C ) gauge theories as a func- 
tion of the number of fermion flavours Nf[l, 2, 3, 4, 5, 6, 7, 8, 9], motivated by the 
pioneering work of Banks and Zaks[10]. For a valuable review, see Ref. [11]. Much of 
the interest stems from the longstanding suggestion that a (nearly) conformal gauge 
theory could play a key role in electroweak symmetry breaking beyond the standard 
model of particle physics[12, 13, 14]. For fermions in the fundamental representation, 
the prediction obtained from a detailed study using the Wilson fermion action[l] dif- 
fers from subsequent studies that use other actions[ll]. According to Ref. [15] it is 
generally believed that confinement exists in the strong coupling limit for any Nf; the 
disagreement with data from the Wilson actionfl] suggests that the unphysical Wilson 
term has a dramatic effect. In the present work, we revisit the strong coupling limit of 
the Wilson action with the goal of providing some additional insight into this issue. 

An SU(N C ) lattice gauge theory with fundamental Wilson fermions is expected to 
have a phase where flavour and parity are spontaneously broken, known as the Aoki 
phase[16, 17]. The expected phase diagram (for example figure 3 of Ref. [17]) is 
sketched in the plane spanned by the two relevant parameters: the gauge coupling 
and the fermion mass. The weak coupling part of the diagram shows "fingers" of 
Aoki phase separated by regions of symmetric phase, but ongoing research by vari- 
ous groups is refining this picture for specific discretizations and lattice improvement 
scenarios[18, 19, 20, 21, 22, 23, 24]. The strong coupling part of the original expected 
phase diagram is simpler: the Aoki phase exists for quark masses of smaller magnitude 
and the symmetric phase exists for quark masses of larger magnitude. In the extreme 
strong coupling limit (g = oo), the phase boundary is predicted to be at a specific 
critical value of the hopping parameter, 



Recall that the hopping parameter k is related to the dimensionless bare fermion mass 
mo as follows: 



In addition, the dimensionless pion mass in the strong coupling limit is predicted to be 
a rather simple function of the hopping parameter [16], 



This expression for the pion mass was derived long ago for quenched configurations (i.e. 
Nf = inside the configurations) [25, 26], and initial tests of its validity for Nf > were 
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presented in some of the direct numerical simulations within Ref. [1]. A central result 
of our work is the observation of systematic deviations from Eq. (3) in exploratory 
simulations with multiple Wilson fermions at ft = 0. 

2. Simulation details 

For this study we use the standard plaquette gauge action and the standard Wilson 
fermion action for each of the Nf degenerate fermions, 
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where ft = 2N C / g 2 with g the gauge coupling, and where U^(x) is the S , [/(A 7 " c )-valued 
link variable. Then the partition function Z is 

, N f 
Z = [dU^x)] \\[di) s (x)][di) f (x)} exp(-S) 
J f=i 

r Hi 

= J \dU,{x)} (detp^ZV)) 2 exp(-,S G ) , (7) 

where D w is the kernel of the single-fermion action S(y = (x)D w (x, y)^(y)- 

Simulations are performed with a standard Hybrid Monte Carlo (HMC). For inde- 
pendent confirmation, simulations for the heavier fermions are repeated using a Poly- 
nomial Hybrid Monte Carlo (PHMC) 1 . In each Molecular Dynamics (MD) evolution, 
the number of steps per trajectory A^md is tuned to beyond 80% acceptance for the 
Metropolis test. The step size At is defined by ATiV MD = 1. 

Our simulations focus primarily on ft = 0, but some studies at ft — 2 will also be 
reported. Lattice sizes include 6 2 x 12 2 , 8 2 x 16 2 , 12 2 x 24 2 and 12 3 x 24. With zero 
temperature simulations in mind, we use periodic boundary conditions and define the 
Euclidean time dimension to be greater than or equal to every spatial dimension. For 
each chosen set of parameters, 50 to 100 configurations were collected after thermal- 
ization, with 4 or 5 trajectories between saved configurations. 

^^Use of our PHMC for the lighter fermions would require a polynomial of remarkably high degree. 
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Our primary observables are the average plaquette (□), the pseudoscalar meson mass 
Itt, and 
fashion: 



m n , and the axial- Ward- Takahashi-identity quark mass m^ WI , all defined in standard 



(□) = /-^^ReTrf/ M (a;)^(a; + ^)f/t( x + ^)^t (x) \ (8) 
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where V denotes the total number of lattice sites, V4 is the symmetric temporal deriva- 
tive operator, and the flavour non-singlet bilinear operators (with flavour index omit- 
ted) are 

P(x,t) = V(x,t)7s^(x,t), (11) 
At(x,t) = ^(x,t) 74 75^(x,t). (12) 



3. SU (2) with many flavours at (3 = 

For a first indication of the difference between quenched and dynamical simulations, 
consider the average plaquette for various numbers of flavours, as shown in Fig. 1. 
An abrupt transition is clearly seen for Nf > 6, with a two-state signal defining the 
transition region in each case. The data at Nf = 4 hint at a transition but are not 
conclusive, though other observables to be discussed below will confirm a Nf = 4 
transition. At Nf = 2 the transition is either suppressed or absent. This behaviour 
implies the existence of a first order transition for sufficiently large Nf, and was already 
reported for the £77(3) case in Ref. [1] and for SU(2) with two flavours of adjoint- 
representation fermions in Refs. [2, 8]. Our results show that the first order transition 
moves to larger values of 1/k as Nf is increased. All of our observed transitions are in 
the range l /& < k < x /i. 

Fig. 2 shows vr? % and m^ WI over a wide range of 1/k. Notice that vr? % and m^ WI 
become insensitive to Nf as 1/k increases. This is no surprise because sufficiently 
heavy quarks have a minimal effect on the vacuum structure of the theory, so meson 
and quark masses for any Nf > approach their quenched (Nf = 0) values as 1/k 
grows. 

A close-up of the small 1/k region is displayed in Fig. 3. The same transition that 
was observed in the elementary plaquette is clearly seen in m\ and m^ WI as well. From 
these data it is clear that the transition exists at least for Nf > 4 and the two-state 
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Figure 1: Average plaquette as a function of 1/k in SU{2) gauge theory with Nf 
flavours at (3 = 0. 

signals imply that it is first order at least for Nf > 6. The transition moves to larger 
1/k as Nf is increased and we speculate that the transition will approach k = y% as 
Nf — > oo, i.e. the gauge theory with infinitely many fermions will become the free 
theory. 

It should be noted that Fig. 1 is obtained for all four of our lattice volumes: no volume 
dependence is observed. This was expected in advance because the average plaquette 
is a short distance quantity. We have similarly confirmed volume independence for 
and m^ WI over our range of volumes, and examples are presented in Fig. 4. In 
particular, the transition shows no evidence of vanishing in the infinite volume limit, 
and it shows no evidence of moving to k = l /i in the infinite volume limit. 

Besides the transition itself, Figs. 3 and 4 show an interesting curvature at values of 
1/k, above the transition, where the slope deviates more and more from the quenched 
slope when approaching the transition point. If someone wanted to extrapolate to 
= using only data above the transition, then the critical hopping parameter so 
obtained, let's name it K^ xt , would be smaller than x /4. Of particular interest is the fact 
that even the data above the transition in Figs. 3 and 4 do not follow Eq. (3) and do 
not arrive at K° xt = x /a despite Eq. (1). 

To ascertain the properties of the transition and of the regions it separates, a range 
of other observables should be discussed. That discussion goes beyond the scope of the 
present work, but brief comments and selected plots can be found in Appendix A. 
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Figure 2: m 2 and m^ WI as functions of 1/k in SU{2) gauge theory with Nf flavours at 
(3 = 0. 

4. SU(3) with many flavours at (3 = 

Qualitatively the situation for SU(3) is similar to that for SU(2), but our quantitative 
results indicate that deviations from Eq. (3) are smaller for SU(3). Perhaps a large- 
N c suppression mechanism is at work. Data for the average plaquette, pseudoscalar 
meson mass, and axial- Ward- Takahashi-identity quark mass are plotted in Figs. 5 and 
6. Notice that m 2 is consistent with Eq. (3) for Nf = 0, but the Nf > data display 
bending that is very similar to the SU(2) case but of smaller magnitude. 

The search for a possible transition is more expensive for SU(3) than for SU(2) 
because we must work at larger k. A transition is found to exist for Nf > 8 but we 
cannot say whether it is abrupt. Recall that Ref. [1] found the transition to occur only 
for Nf > 7; they rely to some extent on counting the conjugate gradient iterations 
required during thermalization to determine the phase for Nf = 6 at k = 1/4. (We 
comment briefly in Appendix B.) Our results 2 do show hints of a transition for Nf = 6. 



2 Thc configurations of Nf = 6 at n = 0.25 in our case are generated from the last configuration of 
Nf = 8 at K = 0.25. We use the HMC algorithm with a time step of At = 0.005, while Ref. [1] uses 
the approximate R algorithm with At = 0.01. 
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5. Discussion 



Simulations on small space-time lattices in the strong coupling limit display the phase 
transition already discussed in Ref. [1], but the present work makes the additional 
observation that the dependence of the pseudoscalar meson mass on the hopping pa- 
rameter deviates from the expected formula, Eq. (3), as seen in Figs. 4 and 6. These 
deviations, which appear as bending in plots of meson (and quark) mass versus hopping 
parameter, are easily seen for N c = 2 and are smaller but still observed for N c = 3. 
The formula in Eq. (3) was derived through a large- N c expansion [16], and the present 
data for N c = 2, 3 suggest a rather rapid approach to Eq. (3) as N c is increased. Eq. (3) 
is also obtained for any N c with Nf = 0. 

Because of the phase transition there is generally no value of the hopping parameter 
at which the pseudoscalar meson becomes massless, but extrapolation from above the 
phase transition (i.e. from smaller hopping parameters) does allow the definition of an 
effective critical hopping parameter, K° xt (iVy). This extrapolated value is found to obey 

< nT(N f - 1) < . . . < < xt (1) < «?* (0) = \ (13) 

as a consequence of the observed deviation away from Eq. (3). 

The authors of Ref. [1] performed a detailed study of SU(3) gauge theory with Nf 
Wilson fermions and from that study they proposed a phase structure with regions 
of confinement and deconfinement for certain values of Nf. They also include a brief 
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Figure 4: Volume independence of m% and m^ WI as a function of 1/k in SU{2) gauge 
theory for Nf = 6 and 12 at (3 = 0. The prediction from Eq. (3) is shown as a solid 
blue curve that reaches = at 1/k — 4. Quadratic fits to lattice data for are 
shown as a pair of solid curves (Nf = 6 is black and iVy = 12 is orange). 

appendix (in the most recent publication in [1]) devoted to the SU(2) theory, where 
k c = 1/4 is assumed and simulations at that hopping parameter are used as the basis for 
conclusions about a confinement / deconfinement transition. The results of Ref . [1] differ 
from those obtained using other actionsfll], and one might ask whether assumptions 
about n c = x /a might play a role. However, Eq. (13) confirms that for any Nf 

K d (N f )<KT(Nf)<\, (14) 

where is the location of the phase transition; this fact is sufficient to leave the con- 
clusions of Ref. [1] intact. Therefore the source of discrepancies between the phase 
structure obtained from the Wilson action versus other actions must be sought else- 
where. 

It is worth noting that phenomena in the strong coupling limit are not typical of re- 
sults at weaker couplings. For example the average plaquette, meson mass, and quark 
mass for the SU (2) theory at (3 — 2 are plotted in Figs. 7 and 8. Instead of an abrupt 
transition that is essentially independent of lattice volume, we now see smoothly con- 
tinuous functions and the meson mass has a significant volume dependence. Detailed 
studies at (3 7^ are left for future work. 
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Figure 5: Average plaquette as a function of 1/ac in S77(3) gauge theory with Nf 
flavours at (3 = 0. 

Appendix A: Additional observables 

A Polyakov loop is the trace of a product of links in a straight line along a lattice axis 
which, via periodic boundary conditions, closes upon itself. Polyakov loops are often 
used as an indicator for confinement. The absolute value of a Polyakov loop is plotted 
in Fig. 9 at (3 = on 6 2 x 12 2 lattices. For each of the two cases Nf = 6 and Nf = 8, 
the graph shows the Polyakov loops in both the shorter (N x = 6) and longer (N t = 12) 
lattice directions. The shorter Polyakov loops are numerically large and dependent on 
the hopping parameter, with a two-state signal at the transition for large enough Nf. 
The longer Polyakov loops are consistent with zero for all values of k. We refrain from 
drawing conclusions about confinement from this brief consideration of Polyakov loops. 

For discussions of chiral symmetry and its possible breaking, natural quantities in- 
clude the chiral condensate (ipip) and the pseudoscalar density (^^ip). The low-lying 
eigenvalues of the Wilson-Dirac operator relate to these observables through the spec- 
tral representation of the quark propagator. For instance, (ip^ip) = x(Hw) an< ^ 

= Yl x ^xjiw) ' wnere ^(Hw) and |A) are an eigenvalue and eigenvector respec- 
tively of the hermitian Wilson-Dirac operator Hw = 'y^Dw. 

The upper panel of Fig. 10 is the lowest eigenvalue, /i, defined by pi = a/A (-£/|^) 
in the eigenvalue equation, H^\\ ) = A |A ), and obtained by the Ritz functional 
method. This behaviour hints at the dependence (^75^) ~ l//i. At least for Nf > 4, 
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Figure 6: m 2 and m^ WI as functions of 1/k in SU(3) gauge theory with Nf flavours at 
(3 = 0. The prediction from Eq. (3) is shown as a solid blue curve reaching m 2 = at 
1/k = 4. 



two phases are found and in both phases the eigenvalues are not particularly small. 

Because Wilson fermions break chiral symmetry explicitly, the chiral condensate is 
not adequate and should be replaced by the subtracted chiral condensate defined by 
(^)subt = 2mf VI /V where TV = (2k) 2 J2 x , t ( P ( x > t ) P ( ' °))- The fermion mass mf* 1 
itself is discussed in the main body of this article, and it is most naturally defined 
only for 1/k values above the transition, so here we focus on the observation of TV. 
The lower panel of Fig. 10 displays Af for SU (2) gauge theory at strong coupling for 
various Nf. Indications of a l/m q divergence in the confinement phase (i.e. 1/k above 
the transition) are visible, as expected. For both observables in Fig. 10, the transition 
is visible when Nf > 4. 



Appendix B: Algorithmic issues 

The HMC algorithm uses a number of trajectories, traditionally distinguished by an 
integer r. Each trajectory uses some number of molecular dynamics steps, ^md, and 
each of those uses some number of Conjugate Gradient (CG) iterations, nf (t) with 
1 < i < Nmd, to attain the imposed precision. The average number of CG iterations 
per MD step is therefore Ncg( t ) — j^— ^^=[ D n ? G ( r )- The left panel of Fig. 11 shows 
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Figure 7: Average plaquette as a function of 1/k in SU{2) gauge theory with Nf 
flavours at (3 = 2 on 8 2 x 16 2 lattices. 

this quantity for sequential trajectories during thermalization for SU(2) with Nf = 6, 
f3 — 0, and k = 0.215 on a 8 2 x 16 2 lattice with precision defined by < 10 -12 . 

For reference, the average plaquette is also plotted for each trajectory. The hot start 
requires more trajectories before reaching thermalization, and some of those trajectories 
needed many CG iterations. One could possibly be tempted to abort such a calculation 
before reaching equilibrium and assume that the increasing Ncg{ t ) would never return 
to a smaller value at larger r. The right panel of Fig. 11 shows Ncg{ t ) f° r £77(3) at 
N f = 6,P = 0, and k = 0.250 on a 6 2 x 12 2 lattice with \\r^\\ < 10" 12 and N MD = 
200. When configuration generation starts from the last configuration of Nf = 8 at 
k = 0.250, the HMC becomes stable at 65% acceptance and Ncg( t ) = O(10 3 ). On 
the other hand, the hot and cold starts are not accepted by the HMC's Metropolis test 
with N C c( r ) = O(10 4 ). In both cases 3 of Fig. 11, the large value of N C g( t ) does not 
imply the existence of a massless mode. 

Acknowledgments 

This work was supported in part by the Natural Sciences and Engineering Research 
Council of Canada, and was made possible by the facilities of the Shared Hierarchical 

3 We plotted the beginning O(10) trajectories in about 300 runs. 



11 



5 
4 
3 
2 
1 


1.5 

1 

0.5 



-0.5 




N f -0 (quenched) 

G-O N r 2 

□ -□ N r =4 
N=12 

N ( =12 (12 2 xl6 2 ) 



-* 



■v" 



5.5 



6.5 

1/K 



7.5 



Figure 8: m 2 and m^ WI as functions of 1/k in SU{2) gauge theory with Nf flavours at 
(3 = 2. Lattices are 8 2 x 16 2 except for one point labelled as 12 2 x 16 2 . 



Academic Research Network (SHARCNET:www.sharcnet.ca). G.K. thanks the Uni- 
versity of Regina for hospitality during her stay there, and we acknowledge use of 
Regina's VXRACK computer cluster. 



References 

[1] Y. Iwasaki, K. Kanaya, S. Sakai and T. Yoshie, Phys. Rev. Lett. 69, 21 (1992) 
and Y. Iwasaki, K. Kanaya, S. Kaya, S. Sakai and T. Yoshie, Phys. Rev. D 54, 
7010 (1996) [arXiv:hep-lat/9605030] and Y. Iwasaki, K. Kanaya, S. Kaya, S. Sakai 
and T. Yoshie, Phys. Rev. D 69, 014507 (2004) [arXiv:hep-lat/0309159]. 

[2] S. Catterall and F. Sannino, Phys. Rev. D 76, 034504 (2007) [arXiv:0705.1664 
[hep-lat]] and S. Catterall, J. Giedt, F. Sannino and J. Schneible, JHEP 0811, 
009 (2008) [arXiv:0807.0792 [hep-lat]]. 

[3] T. Appelquist, G. T. Fleming and E. T. Neil, Phys. Rev. Lett. 100, 171607 (2008) 
[Erratum-ibid. 102, 149902 (2009)] [arXiv:0712.0609 [hep-ph]] and Phys. Rev. D 
79, 076010 (2009) [arXiv: 0901 .3766 [hep-ph]]. 



12 



0.25 



0.2 



0.15 



0.1 



0.05 





1 


1 1 1 


1 1 1 


| 


i | 


i | 

V 


i | i 
N -R (N -6) 


- 


W 













N=6 (N x =6) - 






w 








A 


N=8 (N=12) - 


- 






I 






o 


N=6 (N=12) I 
















- 


















- 
- 








W 






- 
- 


- 

— 






I 








- 

— 










* < 


> * 


¥ 


















V 


















JD 


o 










ft 






1 


1 


1 , 1 


1 


i 


i 1 


1 , - 



3.8 



4.2 



4.4 



4.6 4.i 

1/K 



5.2 



5.4 



5.6 



Figure 9: The dependence of a Polyakov loop (\L\) on Nf flavours is expressed through 



the difference (|L|) 



N, 



\L\)q. This difference is plotted as a function of 1/k in SU(2) 



gauge theory on 6 2 x 12 2 lattices at /5 = with iVj = 6 and 8. 



[4] Y. Shamir, B. Svetitsky and T. DeGrand, Phys. Rev. D 78, 031502 (2008) 
[arXiv:0803.1707 [hep-lat]] and Phys. Rev. D 79, 034501 (2009) [arXiv:0812.1427 
[hep-lat]]. 

[5] A. Deuzeman, M. P. Lombardo and E. Pallante, Phys. Lett. B 670, 41 (2008) 
[arXiv:0804.2905 [hep-lat]] and arXiv:0904.4662 [hep-ph]. 

[6] L. Del Debbio, A. Patella and C. Pica, arXiv:0805.2058 [hep-lat] and 
arXiv:0812.0570 [hep-lat]. 

[7] Z. Fodor, K. Holland, J. Kuti, D. Nogradi and C. Schroeder, arXiv:0809.4888 
[hep-lat] and arXiv:0809.4890 [hep-lat] and arXiv:0907.4562 [hep-lat]. 

[8] A. J. Hietanen, J. Rantaharju, K. Rummukainen and K. Tuominen, JHEP 0905, 
025 (2009) [arXiv:0812.1467 [hep-lat]] and A. J. Hietanen, K. Rummukainen and 
K. Tuominen, arXiv:0904.0864 [hep-lat]. 

[9] A. Hasenfratz, arXiv:0907.0919 [hep-lat]. 

[10] T. Banks and A. Zaks, Nucl. Phys. B 196, 189 (1982). 

[11] G. T. Fleming, PoS LATTICE2008, 021 (2008) [arXiv:0812.2035 [hep-lat]]. 



13 



10"' 

io- 2 
io- 3 

30 
20 
10 



3.8 



— 







A' 



* 



N=0 

0--0 N f= 4 
O-O N =6 

_l I L 



i 1 i r 

— *" N f=° 
N=4 

0-0 N f= 6 
A- -A N=8 



I i L 







A 



4.2 4.4 4.6 4.8 
1/K 



5.2 5.4 



5.6 



Figure 10: SU{2) gauge theory with Nf flavours at (3 = on a 6 2 x 12 2 lattice. Upper 
panel: /i = a/ Ao(-f^y), where Ao(if^) is the lowest eigenvalue of H^, as a function of 
\/k. Lower panel: M = (2k) 2 ^ xt (P(x, t)P(O)) as a function of 1/k. 



[12] B. Holdom, Phys. Rev. D 24, 1441 (1981) and Phys. Lett. B 150, 301 (1985). 

[13] K. Yamawaki, M. Bando and K.-i. Matumoto, Phys. Rev. Lett. 56, 1335 (1986). 

[14] T. Appelquist and L. C. R. Wijewardhana, Phys. Rev. D 35, 774 (1987) and Phys. 
Rev. D 36, 568 (1987). 

[15] T. DeGrand and A. Hasenfratz, arXiv:0906.1976 [hep-lat]. 

[16] S. Aoki, Phys. Rev. D 30, 2653 (1984) and S. Aoki, Phys. Rev. Lett. 57, 3136 
(1986). 

[17] S. Aoki, Prog. Theor. Phys. Suppl. 122, 179 (1996) [arXiv:hep-lat/9509008]. 

[18] F. Farchioni et a/., Eur. Phys. J. C 39, 421 (2005) [arXiv:hep-lat/0406039] . 

[19] C. Michael and C. Urbach [ETM Collaboration], PoS LAT2007 (2007) 122 
[arXiv:0709.4564 [hep-lat]]. 

[20] Ph. Boucaud et al. [ETM collaboration], Comput. Phys. Commun. 179, 695 (2008) 
[arXiv:0803.0224 [hep-lat]]. 



11 




X = trajectory X = trajectory 



Figure 11: Left panel: The average plaquette and Ncg( t ) versus trajectory index r in 
5*7(2) gauge theory on 8 2 x 16 2 lattice at (3 = and k = 0.215 with N f = 6. Right 
panel: N C g(j) versus r in SU(3) gauge theory with Nf = 6 on 6 2 x 12 2 lattice at 
(3 = 0, k = 1/4 and N MD = 200. 

[21] E. M. Ilgenfritz, W. Kerler, M. Muller-Preussker, A. Sternbeck and H. Stuben, 
Phys. Rev. D 69, 074511 (2004) [arXiv:hep-lat/0309057]. 

[22] F. Farchioni et ai, PoS LAT2005, 033 (2006) [arXiv:hep-lat/0509036]. 

[23] V. Azcoiti, G. Di Carlo and A. Vaquero, Phys. Rev. D 79, 014509 (2009) 
[arXiv:0809.2972 [hep-lat]]. 

[24] S. R. Sharpe, Phys. Rev. D 79, 054503 (2009) [arXiv:081 1.0409 [hep-lat]]. 

[25] K. G. Wilson, in New Phenomena in Subnuclear Physics, Proceedings of the 1975 
International School of Subnuclear Physics, Erice, 1975, edited by A. Zichichi 
(Plenum, New York, 1977). 

[26] N. Kawamoto, Nucl. Phys. B 190, 617 (1981). 



15 



